home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cspmv.z / cspmv
Encoding:
Text File  |  1996-03-15  |  3.6 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCSSSSPPPPMMMMVVVV((((3333FFFF))))                                                            CCCCSSSSPPPPMMMMVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CSPMV - perform the matrix-vector operation   y := alpha*A*x + beta*y,
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE CSPMV( UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY )
  13.  
  14.          CHARACTER     UPLO
  15.  
  16.          INTEGER       INCX, INCY, N
  17.  
  18.          COMPLEX       ALPHA, BETA
  19.  
  20.          COMPLEX       AP( * ), X( * ), Y( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      CSPMV  performs the matrix-vector operation
  24.  
  25.      where alpha and beta are scalars, x and y are n element vectors and A is
  26.      an n by n symmetric matrix, supplied in packed form.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      UPLO   - CHARACTER*1
  31.             On entry, UPLO specifies whether the upper or lower triangular
  32.             part of the matrix A is supplied in the packed array AP as
  33.             follows:
  34.  
  35.             UPLO = 'U' or 'u'   The upper triangular part of A is supplied in
  36.             AP.
  37.  
  38.             UPLO = 'L' or 'l'   The lower triangular part of A is supplied in
  39.             AP.
  40.  
  41.             Unchanged on exit.
  42.  
  43.      N      - INTEGER
  44.             On entry, N specifies the order of the matrix A.  N must be at
  45.             least zero.  Unchanged on exit.
  46.  
  47.      ALPHA  - COMPLEX
  48.             On entry, ALPHA specifies the scalar alpha.  Unchanged on exit.
  49.  
  50.      AP     - COMPLEX array, dimension at least
  51.             ( ( N*( N + 1 ) )/2 ).  Before entry, with UPLO = 'U' or 'u', the
  52.             array AP must contain the upper triangular part of the symmetric
  53.             matrix packed sequentially, column by column, so that AP( 1 )
  54.             contains a( 1, 1 ), AP( 2 ) and AP( 3 ) contain a( 1, 2 ) and a(
  55.             2, 2 ) respectively, and so on.  Before entry, with UPLO = 'L' or
  56.             'l', the array AP must contain the lower triangular part of the
  57.             symmetric matrix packed sequentially, column by column, so that
  58.             AP( 1 ) contains a( 1, 1 ), AP( 2 ) and AP( 3 ) contain a( 2, 1 )
  59.             and a( 3, 1 ) respectively, and so on.  Unchanged on exit.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCSSSSPPPPMMMMVVVV((((3333FFFF))))                                                            CCCCSSSSPPPPMMMMVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      X      - COMPLEX array, dimension at least
  75.             ( 1 + ( N - 1 )*abs( INCX ) ).  Before entry, the incremented
  76.             array X must contain the N- element vector x.  Unchanged on exit.
  77.  
  78.      INCX   - INTEGER
  79.             On entry, INCX specifies the increment for the elements of X. INCX
  80.             must not be zero.  Unchanged on exit.
  81.  
  82.      BETA   - COMPLEX
  83.             On entry, BETA specifies the scalar beta. When BETA is supplied as
  84.             zero then Y need not be set on input.  Unchanged on exit.
  85.  
  86.      Y      - COMPLEX array, dimension at least
  87.             ( 1 + ( N - 1 )*abs( INCY ) ).  Before entry, the incremented
  88.             array Y must contain the n element vector y. On exit, Y is
  89.             overwritten by the updated vector y.
  90.  
  91.      INCY   - INTEGER
  92.             On entry, INCY specifies the increment for the elements of Y. INCY
  93.             must not be zero.  Unchanged on exit.
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.